On Metric Preheating 






K. Jcdamzikcl 
Max-Planck-Institut fur Astrophysik, Karl-Schwarzs child- Str. 1, 85740 Garching, Germany 

G. Sig]0 

University of Chicago, Department of Astronomy and Astrophysics, 5640 South Ellis Avenue, Ghicago, IL 60637 

DARG, UMR-8629, CNRS, Observatoire de Pans-Meudon, F-92195 Meudon Gedex, France 

We consider the generation of super-horizon metric fluctuations during an epoch of preheating in 
the presence of a scalar field x quadratically coupled to the inflaton. We find that the requirement 
of efficient broad resonance is concomitant with a severe damping of super-horizon 5x quantum 
fluctuations during inflation. Employing perturbation theory with backreaction included as spatial 
averages to second order in the scalar fields and in the metric, we argue that the usual inflationary 
prediction for metric perturbations on scales relevant for structure formation is not strongly modified. 
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Introduction. In the inflationary paradigm the universe 
went through an early and extended period of accelerated 
expansion a > with rapid growth of the scale factor a. 
Inflation may provide solutions to a variety of cosmologi- 
cal "problems" as well as making generic predictions for 
a scale-invariant Harrison-Zel'dovich spectrum of density 
fluctuations, possibly acting as seeds for structure for- 
mation [yyj. After an inflationary period cosmic energy 
density is dominated by oscillations of the inflaton tf. Ef- 
ficient conversion of this energy to other coupled bosons 
by parametric resonance phenomena may occur during a 
period of preheating immediately following the end of in- 
flation |^,||. Consider, for example, a chaotic inflationary 
scenario with Lagrange density L = ^ip.f^ip''^ — V{(p,x) 
and potential 
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where x is a boson coupled to the inflaton (p. After the 
end of slow roll inflation, i.e. for cp < (127r)~^/^, the 
linearized equation of motion for the comoving Fourier 
mode fc = |k| of the x field is. 



Sxk + SHSjck 
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where the amplitude <I> of the sinusoidal inflaton oscil- 
lations. If — <i>sin(mi), decays as $ = $(io)(a/ao)^'^^^, 
due to the Hubble expansion of the universe which is 
governed by the Friedmann equation 
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for the scale factor a{t), where k^ = 8tt/M^i, with Mpi 
the Planck mass. 

Eq. (Q) is essentially known as the Mathieu equation 
which exhibits parametric resonance, i.e. exponentially 



growing modes for Sxk, due to the time-dependence of 
the X m.ass. One may associate a resonance-parameter 
q{t) — g^$(t)^/4m^ with the system. In an expanding 
universe, parametric resonance is efficient for q > 1 where 
resonance occurs for wavevectors k within broad bands, 
and is rendered inefficient for g ^ 1 Q. In the broad 
resonance regime q > 1, exponential growth of bosonic 
fluctuations is possible for super-horizon modes even in 
the limit fc — > 0. Recently, the important realization has 
been made that growth of coherent, bosonic matter fluc- 
tuations during preheating is accompanied by production 
of metric fluctuations possibly leading to abundant pri- 
mordial black hole formation g . It has also been pointed 
out that a significant modification of the power spectrum 
of density fluctuations relevant for structure formation 
and cosmic microwave background anisotropics is con- 
ceivable |^,q|. 

Multiple scalar fields and metric perturbations. Con- 
sider perturbations around a Robertson- Walker metric in 
longitudinal gauge 



ds^ = (1 + 2(f))dt'^ - a^{t){l - 2(j))dy. 
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where (/> is a gauge-invariant potential, quantifying the 
density perturbation 5p/ p at horizon-crossing in unper- 
turbed Hubble flow. The Fourier-transformed, first-order 
Einstein equations for coupled bosons fi{t) and their 
perturbations 5ipi{t,ii) give [Q 

iH^k + [{k/af + 3i?'] cpk - 

5(pik + SHdipik + [{k/af + V^^^,] 6ipik = 
A^kipi + 2{(pi + 3H(pi)(l)k - ^ V^,^j6(pjk , (6) 

^2 

4>k + H(t>k = -^ 2_/ ^i^^ik ' C^) 



where a subscript Lpi denotes derivative with respect to 
(y9/. Eq. (||) and Eq. (0) may be combined to 

exphcitly showing that (j) is completely determined when 
the evolution of the matter fields is known. 

Whether or not parametric resonance for k —^ in- 
duces significant secondary perturbations in is largely 
a question of initial conditions. The evolution equations 
for the unperturbed, background fields are ifi + SHipi + 
y<Pi — 0. Consider a chaotic inflationary period with in- 
flaton potential as in Eq. (|l|). Assuming that inflation is 
driven by (p (i.e. g^x^ "^ "t-^) ^^ have (p = —{m?l'iH)ip 
by the slow-roll approximation for (p > (/Sq — (127r)^^'^. 
In the limit gip ^ H the evolution of the homogeneous 
X field may be derived in the adiabatic approximation 
X oc \gtp\~^^'^a~^^'^cos J gifdt, such that the amplitude 
of X is decaying during a phase of exponential growth 
of the scalefactor. The adiabatic approximation holds 
if Qq {ip/ipo)~^ <C 2, quite generally when the reso- 
nance parameter at the beginning of preheating satisfies 
go ^ 1- The amplitudes of fluctuations in the x fisld 
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similarly decay as Sxk oc uj~ 
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cos J uj^^dt where 



cj^^ = (k/a)'^ + {g(p)'^ is the oscillation frequency of mode 
k. This can be seen by noting that the R.H.S. of Eq. (0) 
contains factors of x ^^d x whose values are strongly 
suppressed compared to ip and ip due to the inflation- 
ary expansion. When interpreted classically, the decay 
of the amplitude of the Sxk fluctuations corresponds to a 
dilution of particle number densities associated with the 
Sxk field by the expansion of the universe. In contrast, 
super-horizon fluctuations in the inflaton field Sipk are 
not damped, i.e. dpk — constant, and lead to the usual 
prediction of a scale-invariant Harrison-Zel'dovich spec- 
trum of density fluctuations. Hawking radiation in the 
de Sitter phase of an inflationary universe enhances the 
Sipk amplitude for superhorizon modes, which is possi- 
ble since the mass of the inflaton m is somewhat smaller 
than the Hubble constant H. Demanding M^ — gip '^ H, 
in order to have a similar effect for dxk fluctuations, we 
arrive at the condition qq ^ 1/6^' ^. In particular, damp- 
ing seems inevitable for systems in the broad resonance 
regime. 

Consider the amplitudes of vacuum quantum fluctua- 
tions in the bosonic fields on sub horizon scales k/a > H. 
These may be estimated from quantum field theory in 
flat space where the (5(p/'s are replaced by operators 
Sipi{t,x) = /(d3k)/(27r)3[ak(5(?/ke-*'' + H.C.], with at 
the usual annihilation operators and H.C. denoting the 
hermitian conjugate. A typical amplitude may be es- 
timated by the square-root of the expectation value of 
(^j(i,x) in the vacuum. This gives 

(0|5^|(i,x)|0)^ /-^|5^,kP- /#^, (9) 
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(fc/a)^ -I- M^^. Therefore, the vacuum am- 
plitudes in a comoving wavevector range k are given by 
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(10) 



At horizon exit of a mode, i.e. at k/a = H, perturba- 
tions in the inflaton therefore have S(pk\ex — i?/(27r^)^/^ 
since m ^ H. Similarly, for the massive field 5xk we find 
5xk\ex ^ iJ3/V|g(/?|i/V(27r2)i/2. in what follows we wiU 
compute the evolution of the boson fields by the classical 
equations which is a good approximation in the limit of 
large occupation numbers, in particular, when the am- 
plitudes are much larger than their vacuum expectation 
values. 

Let us evaluate the wavevectors for modes which con- 
tribute to the formation of large-scale structure. Let us 
define, Tq = (mipo)^'^, which is the approximate reheat 
temperature, Thh — To, if reheating would be instanta- 
neous. For what follows it is sufficient to approximate 
reheating to be instantaneous, such that a evolves as a ~ 
(t/toy/^ for t :$> to [setting a{to) = 1]. Then the proper 
horizon distance evolves as r^(r) - r^^ (To)(To/T)2 
with r^(To) ~ Mpi/Tg, the horizon distance immedi- 
ately after inflation. A mode with comoving wave vector 
k ~ {Tre/To)/rPj{To) ~ {m/3){Tre/To) re-enters into the 
horizon at T^e <C Tq. This mode left the horizon dur- 
ing inflation at time tex when k/aex — H(p>ex)- This 
allows one to derive the scale factor at horizon exit, 
aex — (Tre/To){ipex/po)~^- For k/a < gip the ampli- 
tude of the matter field at the beginning of preheating 

Sxk{to) = 6x{tex)alx{Pex/poy^'^ IS givcu by 



5xk{to) ^ O.lmgg ^' {Tre/Tof^'^{pex/poy 



(11) 



For small k and low Tre ''^ leV this is very small com- 
pared to a typical inflaton perturbation Spk ^ "m- 

Non-perturbative decay of the inflaton condensate into 
bosonic matter fields via parametric resonance occurs 
when the condition of adiabaticity for the x field is vi- 
olated, i.e. when iu-^^. ^ w^ . In what follows, we fol- 
low closely the analytic treatment by Kofman, Linde, 
and Starobinsky Q (KLS97) for the stochastic growth 
of bosonic particle number in the broad resonance pa- 
rameter regime, g ^ 1, and in an expanding universe. 
The largest bosonic wavevector kmax which receives para- 
metric amplification is given by kmax/o- — (.g?n<I>)^"/2. 
Due to phase space arguments, this is also the typical 
wavevector k^ in which most of the inflaton energy is 
"dumped" . Let us define 5Xk as the amplitude of the 
boson fields, i.e. 5Xk = 5xk at moments when 5xk — 0. 
KLS97 were able to derive that the comoving occupa- 
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tion number n^ evolves as, 2nk — ^Sxk \^^i 
exp(277i J di/ifc(t)), where the integral may be replaced 
by 2niiit, with an effective value ^ ~ 0.1 — 0.18. Re- 
placing the Floquet index by an effective index accounts 
for the stochasticity of the resonance, due to wavevectors 



shifting in and out of resonance bands with the expan- 
sion of proper wavelength and the time-dependency of 
the resonance parameter q{t). Note that there is still 
some dependency of // on fc, which implies that the de- 
cay of the inflaton energy is dominated by decay into 
modes with fc in a narrow neighborhood of the maximum 
unstable mode. These considerations allowed KLS97 to 
estimate the number of inflaton oscillations N — mt/{2-K) 
between t^ ^ t ^ ti, where ti is the time when backre- 
action effects of the produced matter particles on the 
inflaton become significant, i.e. when m^ ~ 5^X^„^(ii), 
where in the following for any function /(x) we define 
/™, = (/2(x))i/2. They found, 
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During this time the scale- factor evolves as a{t) = 
{i/ioY^^ — (ttA^)^/^, such that at ti the resonance pa- 
rameter is qi = g'^^'^{ti)/{Am?) ~ qo/(7riVi)^- Subse- 
quently, if qi > 1/4, the inflaton oscillation frequency 
increases to m? ^g ~ g^X^^^, and Hubble expansion ef- 
fects are negligible, i.e. m^^cS ^ H. The second stage of 
preheating, t > ti, is well approximated by parametric 
resonance according to the Mathieu equation with res- 
onance parameter qcs — .9^^^/(4(?^^rms)- Broad reso- 
nance is terminated at some time ^2 , when q^s decreased 
to ~ 1/4. At this time $^(^2) - X^^^{t2). During the 
second stage of preheating, the inflaton typically makes 
only a few oscillations, A^2 — (l/47r/i) ln4g]^ . Combin- 
ing adiabatic evolution with resonant ampliflcation for 
wavenumber k between tg and ^2 yields 



SXk{to) 
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exp [27r^(iVi + N2)] 
a3/2(t2) 



(13) 



The contribution of x fluctuations to the large-scale 
metric fluctuations at ^2, if estimated via Eqs. (^ and 
( |l3| ) and employing the severely damped background x, 
would be utterly negligible. Nevertheless, the problem 
is intrinsically non-linear since Sxk S> X- If instead, we 
were to include the second order contribution such as ~ 

(.gvs)^ /d^k'(5^/k"^^/(k-k')' 01^ t^'3 R.H.S. of Eq. (||) we 
would flnd significant metric perturbations when back- 
reaction becomes important, even in the limit A: — > 0, 
sourced by the sub-horizon fluctuations in Sx- This is 
not surprising though, at onset of backreaction there is 
approximate equipartition of energy between the inflaton 
and the coupled boson. Since the energy in the Sx fluctu- 
ations is not taken into account in Eq. (p), one expands 
effectively around an unphysical Hubble constant. We 
may rather attempt to estimate 4>k~,Q{t2) by associating 
an effectively homogeneous background of the x field on 
large scales by taking the root mean square Xj-ms over 
small-scale fluctuations. This would yield (/)fe^olx(^2) ^ 
52$2X™,(5Xfc/52$2^2^^ _ SXk/Xrms (^2), for the con- 
tribution of X fluctuations to the metric perturbations. 



Using Eq. (y) with T^ ~ leV, Tq ~ 4 x 10'^ Mpi, 
and ipex — SMpi we find 6Xk{to) ^ 10^^^. Such a 
small initial amplitude is only amplified by modest fac- 
tors - 10-* - 10^ estimated via Eq. (|l|) for fi ~ 0.13 
and g ^ 3 X 10^^ — 10~^.Nevertheless, we will see below 
that this result represents a gross underestimate. 

Second order contributions to the energy density from 
scalar field and metric fluctuations on scales exceed- 
ing the typical fluctuation scale may be accounted for 
by the effective energy-momentum tensor formalism of 
Ref. pj. This corresponds to substituting the poten- 
tial V = V + ^ J2ij ^•pi'Pj^fi^fJ s-nd adding the term 
[dtp] + a~^(V(5(p/)^]/2 -I- 2V^j4)5ipi for each scalar field 
under the parentheses in Eq. (|^), and using the modified 
equation 



{Cpi + iH^i) (1 + 4(/.2) + V^, ~ 2<j)5(pi - Ac^S^i 
-QHipSipi + 4:ipj(j)<j> - 2a'''^(j)V^5ipj = 



(14) 



for the homogeneous fields. In these expressions second 
order terms are to be understood as spatial averages. 
Whereas for small scales ~ k^\^ such a scheme may 
be questionable, it should be appropriate for the super- 
horizon scales relevant for structure formation. 
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FIG. 1. Results of a numerical simulation of metric pre- 
heating with two scalar fields and parameters as discussed in 
the text. The top panel shows, from top to bottom, the total 
I^Xrmsl, the homogeneous field |x|, and the large scale fluc- 
tuations I^Xfc^ol- The curve for |5xrms| starts at its vacuum 
expectation value at the beginning of preheating. The bot- 
tom panel shows the total |(?!)rms| (dashed lines) and the large 
scale metric potential |0fc_,o| (solid lines) as functions of time. 
The results obtained from integrating Eq. (a) (upper set of 
curves) and Eq. iXm (lower set of curves) diverge once backre- 
action sets in, reflecting the uncertainties within the current 
numerical approach. 

Numerical Model. We have numerically integrated 



Eqs. (||) and and (Q) along with Eq. (|lj) for the homo- 
geneous fields and the Friedman equation (||) with back- 
reaction terms for the two field case and the potential 
Eq. (||). We started integration deep in the inflationary 
regime where the large scale structure modes relevant to- 
day were still within the Hubble horizon. We used the 
initial conditions </? — 4:AIpi, x = "^^/g and Eq. ( pX| ) for 
the matter fields, as well as 5(fik ^ lOikSfik for fluctua- 
tion derivatives (the results are qualitatively insensitive 
to the exact values). The initial value of (pk was evalu- 
ated from Eq. (ph. A grid of 100 equidistant logarithmic 
momentum modes was implemented between k = m/2^/'^ 
and k = 1.5[gTO$(io)]^^^ = 3TO(go/4)^^'* which covers the 
dominant resonant modes around the Hubble horizon. In 
addition, one mode with comoving scale of ^ 100 Mpc 
today, corresponding to T^e — 0.5 eV in Eq.(|ll|), was fol- 
lowed. 

In Fig. |l| we show results of simulations with m = 
10~^Mpi and qo = lO''. Note that Fig. includes re- 
sults from two simulations (a) as outlined above and (b) 
as above but by integrating Eq. (||) instead of Eq. (0). 
In the linear limit (i.e. mt/{2TT) ^ 25) the results of 
both simulations coincide, providing a consistency check 
of our numerical routine. During the inflationary phase 
one may observe the standard growth of metric fiuctu- 
ations due to the increase of (1 -t- it;), where w is the 
ratio of pressure and density. For mt/(2Tr) > 5 para- 
metric resonance develops and for mt/ (27r) > 25 back- 
reaction becomes important. The saturation of X^ms at 
~ m/g = Qq ^{to)/2 ~ 3 X 10^^ due to backreac- 
tion is clearly visible in Fig. |l|. Note that due to the 
choice of go and m, Xrmsit2) — Xrms{ti) since qi ~ 1/4, 
thus the second stage of preheating is very short. When 
backreaction becomes important, results obtained from 
integrating Eq. (0) and Eq. (^, respectively, which are 
equivalent in the linear regime, start to differ. In par- 
ticular, whereas Eq. (ph predicts a second rise in small 
scale metric fluctuations at ^2, and oscillatory behavior, 
integration of Eq. (R) shows a smooth transition. These 
differences are likely due to the replacement of small scale 
fluctuations by numerical phase space sums which may 
introduce some coherence between the perturbed quan- 
tities which in reality, in particular for large scales, does 
not exist. Within this approximative scheme, c^fe^olxC^z) 
is still < 0fc^o|ip(^2), but much larger than our "naive" 
analytical estimate. The coupling to small-scale metric 
(f> and boson perturbations 6x in Eq. (|lj) tends to lift 
the homogeneous x field to ~ (pSX which in turn in- 
duces a source term for Sxk->o in Eq. (^. These source 
terms induce a growth of x and Sxk^o, with growth 
due to parametric resonance subdominant. In contrast, 
the employed model approximates growth of Sxk~k^am > 
as well as backreaction on ip, by the usual parametric 
resonance theory Q, i.e. essentially unmodified by the 
existence of metric fluctuations. Our conclusions, i.e. 



(f)k->o\x ^ 4>k^o\ip and (j)k^o ^ 1, are not likely modi- 
fied by a subsequent epoch of rescattering and thermal- 
ization, but may be if </> grows further by gravitational 
instability. Nevertheless, for ~ 1 efficient formation of 
primordial black holes seems inevitable and the notion 
of an effectively homogeneous and isotropic universe on 
large scales, described by a FRW metric, is called into 
question. 

Summary. We have investigated preheating after an 
inflationary epoch for a simple chaotic inflationary model 
with the inflaton coupled to another bosonic particle. 
Our study focuses on the possibility of generation of 
secondary metric perturbations on scales relevant for 
structure formation. Within first-order perturbation the- 
ory, but with backreaction included to second order in 
the scalar fields and in the metric in an effective way, 
we found that metric fiuctuations on scales relevant for 
structure formation stay well within the linear regime 
during preheating. This is mainly due to the fact that at 
the beginning of preheating the homogeneous component 
of the scalar field x, and its super-horizon quantum fluc- 
tuations Sx, which couple to the inflaton, are severely 
suppressed compared to the inflaton B. Our conclu- 
sions are similar to those of Refs. |^,|l^ who investigated 
the case of general-relativistic parametric resonance of 
perturbations around the coherent oscillations of a mas- 
sive, but otherwise uncoupled inflaton. However, we cau- 
tion that our study does not fully address the possibil- 
ity of efficient mode-mode coupling between sub-horizon 
and super-horizon metric fluctuations. This issue may, 
in principle, be resolved by numerical simulation as in 
Ref. 0. 
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